The Canted Spiral: An Exact Ground State of XXZ Zigzag Ladders 
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We derive the exact ground states for a one dimensional family of S = 1/2 XXZ Hamiltonians on the 
zigzag ladder. These states exhibit true long range spiral order that spontaneously breaks the U(l) invariance of 
the Hamiltonian. Besides breaking a continuous symmetry in d — 1, this spiral ordering has a ferromagnetic 
component along the symmetry axis that can take any value between zero and full saturation. In this sense, our 
canted spiral solutions are a generalization of the SU(2) Heisenberg ferromagnet to non-zero ordering wave- 
vectors of the transverse spin components. We extend this result to the d — 2 anisotropic triangular lattice. 

PACS numbers: 72.80.Sk, 74.25.Ha, 73.22.Gk 



The search for chiral phases in one dimensional frus- 
trated magnets has been very active during the last ten years 
H] |2 13 |U |5] |6). This interest was triggered by the pre- 
diction of a ground state with non-zero vector spin chirality, 
(Sj x Sj+i) + 0, for the J x , J 2 XXZ chain with |J X | < J 2 
(Ji and J 2 are the nearest and next-nearest exchange inter- 
actions). The Jx, J 2 chain is equivalent to a zigzag ladder 
(see Fig{T^) and it has also attracted a lot of interest during 
the last decades QE). One of the main reasons for the con- 
tinuiung fascination generated by this model is in the inter- 
play between geometric frustration and strong quantum fluc- 
tuations that leads to a rich and exotic variety of physical phe- 
nomena. 

The vector chirality has to be distinguished from the scalar 
chirality, (Sj-_i • Sj x Sj+i), which breaks different discrete 
symmetries and is associated with different physical quantities 
[|9l . As it is mentioned in Ref.flU, classical states with sponta- 
neously broken chirality only exist together with helical long 
range order. While the vector chirality distinguishes left and 
right spirals, the sign of the scalar chirality distinguishes be- 
tween positive and negative canting angles for a given spiral 
orientation. The helical order breaks the continuous symme- 
try of global spin rotations along the z-axis. Consequently, 
the existence of long range helical order is in most cases pre- 
cluded by zero point fluctuations of d = 1 quantum systems 
(Mermin- Wagner theorem ifTOl ). On the other hand, chiral 
orderings are allowed because they only break discrete sym- 
metries. It is for this reason that chiral orders in quantum spin 
systems can be thought as remnants of the helical order in 
classical systems. This has been one of the main motivations 
for finding chiral orders in quantum spin Hamiltonians whose 
ground state exhibits helical order in the S — > 00 limit. 

The zigzag XXZ ladder is one of the simplest spin mod- 
els whose classical counterpart has a helical ground state in a 
certain region of exchange parameters. The original proposal 
of a chiral ground state for the quantum version of this model 
was based on a mean field treatment of the bosonized Hamil- 
tonian [1|. A similar approach was used by Kolezhuk and 
Vekua |4| to obtain a field induced chiral state in the Heisen- 
berg (XXX) zigzag ladder. These mean field approaches were 
later validated by numerical simulations [5, 6|, but there are 
some regions of the quantum phase diagram were the situation 



is still unclear (6). In general, the numerical detection of these 
phases is very challenging due to the lack of commensuration 
between the dominant wave-vector of the spin-spin correlator 
and the finite size chains that can be solved numerically ifTTI . 
It is therefore desirable to find exact analytical solutions for 
the chiral ground states and compute the most relevant corre- 
lation functions without making any approximation. 

In this Letter we consider a one dimensional family of 
S=l/2 XXZ zigzag ladders for which we derive an exact 
ground state subspace with a very unusual property: it is the 
semi-classical version of a canted spiral phase (see FigJTJ?). 
In particular, this implies that the ground state has true long 
range helical order, i.e., it breaks spontaneously the contin- 
uous U(l) symmetry of global spin rotations along the z- 
axis. Moreover, the ground state is still a canted spiral on 
finite size ladders for a discrete set of ratios, J1/J2, between 
the two exchange constants. The most curious of aspect of 
these solutions is that in contrast to the usual example of the 
SU(2) Heisenberg ferromagnet (the SU(2) symmetry is also 
spontaneously broken at T = 0), the order parameter of the 
canted spiral does not commute with the Hamiltonian. We 
will see that this observation is related to the emergence of 
an SU(2) symmetry in the ground state subspace lfT2l . The 
nonzero scalar and vector chiralites are an epiphenomenon of 
the helical solution. They are also computed in an exact way 
by exploiting the semi-classical or unentangled nature of the 
solution. We also provide an exact calculation of a branch 
of gapless single-magnon excitations with a quadratic low- 
energy dispersion around the spiral wave-vectors ±Q. Fi- 
nally, we extend our exact canted spiral solution to the case 
of an XXZ Hamiltonian defined on a two dimensional version 
of the zigzag ladder (see Fig[T};). This is indeed the magnetic 
lattice of the frustrated quantum magnet Cs 2 CuCl4 fl3l . 

We start by considering an XXZ model on the zigzag ladder 
with L sites (0 < j < L — 1) depicted in Fig{T] 

where v — 1,2, < j < L — 1, L = (periodic boundary 
conditions), and Sj = Sj ± iSj. We will assume that the 
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FIG. 1: a) Zigzag ladder. The dotted lines indicate that we are using 
PBC's. The arrows show the circulation of the spin and the electrical 
orbital currents that result from the nonzero vector and scalar chi- 
ralities of the canted spiral solution, b) Canted spiral solution. The 
arrow at site j corresponds to (^e.^o l^i l^.vo)- c) Extension of the 
zigzag ladder to d = 2. ei and &2 are primitive vectors. r+ and r_ 
denote the triangular plaquettes that are above and below the site r. 



We claim that the set of linearly independent states 



(6) 



are exact ground states of H for < p < L. The normal- 
ization prefactor, (^f p \^f p ) = 1, is obtained by noticing that 
L p / 2 (Sq) p \%) is a linear combination of (^) states with coef- 
ficients that have the same absolute value equal to pi. 

There are different ways of demonstrating our claim. Here 
we will use an algebraic procedure that unveils an emergent 
SU(2) symmetry of H. In first place, we will demonstrate 
that the states | are degenerate eigenstates of H . For this 
purpose it is convenient to express H in momentum space: 



i 



(7) 



where 



J„ 



JiAi cosq + J 2 A 2 cos2g, 
Ji cos q + J 2 cos 2q. 



(8) 



ratio of exchange constants satisfies the equation 



cos Q 



~4J 2 



(2) 



where Q = 2irn/L (0 < n < L — 1) is a wave-vector in the 
Brillouin zone (BZ) of the ladder. We will also assume J 2 > 
0, so Eq.([2| defines a set of L/2 + 1 different ratios J1/J2 for 
even L (the ratio is the same for ±Q). In the thermodynamic 
limit, L — > 00, Q can take any arbitrary value between and 
7T and Eq.([2| always has a solution for | J\\ < 4J 2 . 

We are particularly interested in the one dimensional family 
of Hamiltonians defined by — A V J V with 



From this expression we derive the following commutation 
relation by using Eqs.Q: 



with 



VL l = l,L 



a l= J2 JvSm"Q(S? +u -S?_ u ). 

v=l,2 



From Eq.([T]i, it is evident that 

#10) 



0. 



(9) 



(10) 



(11) 



Ai — cos Q = 
A 2 = cos 2Q = 



4J 2 



8Jf 



Ground States .The next step is to provide a subspace of exact 
ground states for this Hamiltonian family. For this purpose we 
will define the vacuum, |0) = as the fully polarized 

state with all the spins down: S||0) = -(1/2) |0) V < 
j < L — 1. In addition, we introduce the spin operators in 
momentum space: 



s z q = 



1 



L ^ 

j 



1 



(4) 



and S q = (S+y. These operators obey the following com- 
mutation relations: 



L 



sr. 



1 



k+q- 



(5) 



In addition, our definition of the vacuum state, 5| 
-(1/2)|0) V < j < L - I, implies that (5| - Sf 



(3) V < j,l < L- 1, anda ( | 



0. Consequently, 



[S+H}\®) = 0. 



(12) 



Finally, by using Eqs.Q and |9]) it is easy to verify that 

[S+, [5+ fl]]|0)=O. (13) 

The combination of Eqs.([T2]) and ( [T3] > leads to a very impor- 
tant result: 



[S+ fl](S+)"|0) = O. 



(14) 



This result implies that Sq is an infinitesimal generator of a 
continuous symmetry of H when H is restricted to the sub- 
space Qq generated by the states \^ v ) with < p < L: 



[S+,Pg Q HPg Q }=0, 



(15) 
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where Pg Q is the projector on the subspace Qq. Moreover, 
Eqs.jTTJ and jl4| imply that H (5+) p |0) = or 



HVS/r. 



V <p< L. 



(16) 



This concludes our first step. The second step is to demon- 
strate that the eigenstates \^f p ) are ground states of H. Since 
the corresponding eigenvalues are equal to zero, we just need 
to demonstrate that H is semi-positive definite. This can be 
done by rewriting the Hamiltonian in the following way: 



H 



(17) 



with Hj = h 2 (j - l,j + 1) + hi(j - l,j) + hi{j,j + 1) 



h 2 (l,n) = j*(S?S? l -\) + ^(S+S-+SrS+) 



MU) = J ^(SfS z n -^) + J ^(S+S-+SfS+) (18) 



For Ai and A2 given by Eqs.Q, we obtain that 

Hj = ( J 2 + J!/8J 2 )Pj 



(19) 



where Pj is a projector on the two-dimensional subspace 
£(j — li iii + 1) generated by the states: 

& T ) = 7 (2 cos QS+_ V S+ +1 S+S+ +1 S+^S+Mj, 
= "f(2cosQS+-S+_ 1 -S+ +1 )m^ (20) 



with 7= 1/V4 + 2 cos 2Q and |0>j = ||>j®|| 
Given that J2 > 0, the combination of Eqs.(17i and (19i 



implies that H is semi-positive definite. This concludes the 
demonstration of our main claim: the states \^ p ) generate a 
ground state subspace Qq of H. For Q ^ 0, 7r, there are two 
ground state subspaces, Qq and Qq (Q — —Q), that corre- 
spond to right and left spiral solutions as we will see below. 
In addition to these solutions, there are other linearly inde- 
pendent ground states in the S z = subspace that will be 
presented elsewhere fl4l . 

Order Parameter. The exact ground subspace Qq contains 
a spiral order parameter that breaks spontaneously a continu- 
ous symmetry of H: this is the U(l) symmetry of global spin 
rotations around the z-axis. To see this we just need to choose 
an appropriate linear combination of the |^ p ) ground states: 



I* 



6, Vol 



(21) 



To verify that \^e,ip ) £ Qq we just need to notice that Qq is 
an invariant subspace of Sq and Sq, something that will be- 
come evident when we discuss the emergent SU(2) symmetry 
of H. It is easy to verify that \^e, vo ) = T Qt(po Rg\^ L ) with 



Re 



^(s + -s-) ; T Qm = e^M+nW . (22) 



Re is the global spin rotation by an angle 9 along the y-axis, 
while Tq iVo is the "twist" operator with momentum Q plus 



a global rotation by an angle (po along the z-axis. In other 
words, Xq iVo rotates the spin j by an angle ipj = Qj + cpo 
along the z-axis. This implies that \^e,ip ) is a canted spiral 
solution in wich the spin j is fully polarized along the direc- 
tion Uj = (— sin 8 cos tpj , sin sin ipj , cos 8) (see FigjTJ?): 



S, • 



= 2^0) 



This implies that 1^0, «„) is a direct product state: 



I* 



(23) 



(24) 



j=0,L-l 



where 



l^ob = e l ^cos0/2|T),+e-' i ^sin0/2||) J , (25) 

Besides breaking the continuous U(l) symmetry of H, this 
solution breaks two discrete symmetries: spatial inversion, X, 
and the product of time reversal times a 7r-rotation along the 
z-axis Te ilTs/Zs S. The more physical implication is that the 
canted spiral carries non-zero vector, Kj = Sj x Sj+i ■ z, and 
scalar, Xj — S 3 -_i ■ Sj x S-j+i, spin chiralities: 

(^e, V0 \^j\^e, V0 ) = ^sinQsin 2 6», 

(*e, w \Xi\*e,<P ) = cos0sin 2 0sin 3 (Q/2)cos(Q/2).(26) 

These identities are easily obtained by exploiting the . direct 
product form of \^e, Vo ) [see Eq.(24i]: (Sj x Si) = (Sj) x 
(Si). The coexistence of both chiralities implies that there 



are spin and electric currents circulating around each of the 
triangular plaquettes (see Fig{T^) I91 IT51 . 

Emergent SU(2) symmetry. A spontaneously broken contin- 
uous symmetry is uncommon for d = 1 systems because it is 
usually prohibited by the Mermin- Wagner theorem [ 10 1. The 
simplest counter-example is the SU(2) ferromagnet (Q = 0). 
In that case the continuous SU(2) symmetry is spontaneously 



broken because the order parameter, v L(S§, Sq,S§) , coin- 
cides with the infinitesimal generators of the SU(2) symme- 
try group, i.e., it commutes with the Hamiltonian. Note that 
S* = (S+ + S~)/2 and S«f = i(S~ - S+)/2. The sit- 
uation is less clear for the canted spiral under consideration 
(Q 7^ 0) because the order parameter Pq = (Pq, Pq, Pq) = 

\^L(Sq,Sq,Sq) does not commute with H. However, we 
will show below that the components of the spiral order pa- 
rameter generate an SU(2) group that is an emergent symme- 
try of #Q2. 

The fact that the components of the canted spiral order pa- 
rameterPg = (P^P^P^) = VL(S^,S y Q ,S^) are ele- 
ments of an su(2) algebra results from the commutation rela- 
tions Q: 



[pirn 



(27) 



where e w „ are the components of the Levi-Civita tensor. Ac- 



cording to Eq.( 15 1, Pq commutes with Vg Q HVg Q and this 
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implies that the the SU(2) group generated by Pq is an emer- 
gent symmetry of H [12|. Moreover, the ground state sub- 
space Qq is an irreducible representation of this SU(2) group 
with an eigenvalue of the Casimir operator P • P|^ p ) = 
L(L/4 + l/2)\^ p ) VO < p < L. This eigenvalue coincides 
with the square of the total spin for the uniform Q = case. 

Low Energy Excitations. The U(l) and the translational in- 
variance of H imply that the single and two-magnon excita- 
tions on top of the fully polarized ground states, |0) and I'I'l), 
can also be computed in an exact way. In particular, the wave 
function for the single magnon with wave-vector q on top of 
the fully polarized state |0) (the other one is obtained by ap- 
plying the time reversal symmetry) is 5+|0), and the corre- 
sponding energy eigenvalue is 



J v (cos vq — cos uQ) . 



(28) 



v=1.2 



Two dimensions. Eq.([T9]) shows that H is a sum of projec- 
tors Pj over all the triangular plaquettes of the zigzag ladder. 
This structure suggest a natural extension of H to the d = 2 
case of parallel chains coupled by zigzag bonds (FigjTJ:): 

H 2d = H% + H 2 r d = (J 2 + J 2 /8J 2 ) ^(P r+ + P r- )> 

r r 

(29) 

where 

H 2 r t = h 2 (r±e 1 ,r±e 2 ) + h 1 (r±e 1 .r) + h 1 (r,r±e 2 ) (30) 

and P r ± is the projector on the subspace £(r±ei,r,r±e2) 
[see Eq.(20i]. The sites r belong to the L x L/2 lattice 
depicted in Fig[l]:. Note that H 2d is an XXZ Hamiltonian 
with exchange interactions Jf , Ji along the diagonal bonds 
and 2J|, 2J 2 along the horizontal bonds (now each horizon- 
tal bond is common to two triangular plaquettes). Again, 
the boundary conditions are periodic. The three sites (r ± 
ei,r,r ± e 2 ) are the corners of the triangles that are above 
(+) and below (— ) the site r (see FigfTJ;). 
The d = 2 version of the canted spiral state is 



l*&a> 



\i>e,<p )r, 



(31) 



where \ipg u> a ) T is given by Eq(25 1 with ipj replaced by tp r 
ip + Q • r. Q is a wave-vector of the BZ that satisfies: 



Q e 2 = Q ei = 



(32) 



Again Eq.^ holds for Q = 2m/ L. To prove that 1*^) 
is a ground state of H 2d we simply use the right hand side of 
Eq.(j29Jl and note that \^e, Vo )r± ei ® \4>e, Vo )r ® |V'e, Vo )r±e 2 
belongs to E - L (r±ei, r, r ±e2). This argument also provides 
an alternative way of proving that \^g jV0 ) is the ground state 
in the d = 1 case of H . 

Conclusions. In summary, we have demonstrated that 
the canted spiral solution is the ground state of a family of 
5=1/2 XXZ Hamiltonians defined on the zigzag ladder and 



parametrized by the wave-vector Q of the spiral. The sponta- 
neous breaking of a continuous symmetry is very unusual for 
d = 1 quantum Hamitlonians when the order parameter does 
not commute with H. We have shown that the case under con- 
sideration has an emergent SU(2) symmetry generated by the 
components of the canted spiral order parameter. 

The family of exact canted spiral solutions includes the 
usual SU(2) ferromagnet in the Q = limit. Like in the fer- 
romagnetic case, the Q ^ solutions are product states (un- 
entangled), meaning that inter-site quantum fluctuations that 
typically remove the long range spiral order [ 10 1 are absent 
in this case. Moreover, the exact low-energy spectrum of sin- 
gle magnon excitations exhibits a quadratic dispersion around 
the ordering wave-vectors ±Q. In this sense, we can say that 
our canted spiral ground states are natural extensions of the 
ferromagnetic solution of isotropic Heisenberg models. Note 
also that our family of quasi-exactly solvable Hamiltonians 
includes cases that are arbitrarily close to the isotropic limit 
for < \Q\ <C 7r (Ji < —4J 2 ), i.e., the exchange uniaxial 
anisotropy is a small perturbation. This regime corresponds 
to a ferromagnetic nearest-neighbor interaction J x < ( J 2 
is always antiferromagnetic) that has also been considered in 
previous works lfl6l [T71 [Til [T9l . 

We have also shown how to extend our exact solutions to 
the two dimensional case of XXZ Hamiltonians defined on 
a lattice that comprises horizontal chains coupled by zig-zag 
bonds (see FigjTJ;). This is just a spatially anisotropic trian- 
gular lattice that appears in real quantum magnets such as 
CS2CUCI4 1131 . By showing that H can be written as a sum 



of projectors Pj over triangular plaquettes [see Eq.(19i], we 



are giving a prescription for constructing other higher dimen- 
sional Hamiltonians with exact canted spiral ground states. 

This work was carried out under the auspices of the NNSA 
of the U.S. DOE at LANL under Contract No. DE-AC52- 
06NA25396. 
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